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Abstract 

s 

, We use Faddeev's decomposition to solve the shell-model problem for three 

^ , nucleons. The dependence on harmonic-oscillator excitations allowed in the 

Q I model space, up to 32?ir2 in the present calculations, and on the harmonic- 

■ oscillator frequency is studied. Effective interactions derived from Nijmegen 

• II and Reid93 potentials are used in the calculations. The binding energies 

obtained are close to those calculated by other methods. The structure of the 
Faddeev equations is discussed and a simple formula for matrix elements of 
I the permutation operators in a harmonic-oscillator basis is given. The Pauli 

' principle is properly treated in the calculations. 

o 

l> 
5 , 

I. INTRODUCTION 

V ' 

^ . Many different methods have been used to solve the three-nucleon problem in the past. 

P5 ! The most viable approach appears to be the Faddeev method [0]. It has been successfully 
applied to solve the three-nucleon bound-state problem for various nucleon-nucleon poten- 
tials p|,^,§J^. The most complex calculations of this kind include up to 34 channels, when 



X 



;_j ■ all the j < 4 waves are taken into account. The precision achieved in these calculations is 



better than 1% gi. 

On the other hand, when studying the properties of more complex nuclei one typically 
resorts to the shell model. In that approach, the harmonic-oscillator basis is used in a 
truncated model space. Instead of the free nucleon-nucleon potential, one uses effective 
interactions inside the model space. Examples of such calculations are the large-basis no- 
core shell-model calculations that have recently been performed . In these calculations 
all nucleons are active, which simplifies the effective interaction as no hole states are present. 
The effective interaction is determined for a system of two nucleons in a harmonic-oscillator 
well interacting by the nucleon-nucleon potential and is subsequently used in the many- 
particle calculations. 
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Republic, 250 68 Rez near Prague, Czech Republic. 
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In the present paper we combine the apphcation of the shell-model approach to the 
three-nucleon system with the Faddeev method. We use Faddeev's decomposition for the 
basis states and remove the center-of-mass term. This leads to a significant simplification 
of the problem and allows us to extend considerably the model space. We can study the 
convergence properties of the results with the increasing model space. If convergence is 
achieved, then this approach leads to the exact solution of the three-nucleon problem and is, 
thus, complementary to the traditional calculations, based on either the differential equation 
solutions in the configuration space or integral equation solitions in the momentum space 
treatments. 

In addition to the attempt of solving the three-nucleon problem exactly, the present 
method serves primarily as a test of the shell-model approach. In particular, it allows 
us to test effective interactions used in standard shell-model applications. The present 
approach has several advantages. First, any number of partial waves can be included. 
Second, the calculation is simplified by using a compact formula for the matrix elements of 
the permutation operators in the harmonic-oscillator basis. Also, because of the way we do 
the model-space truncation, we keep equivalence of the Faddeev and Schrodinger equations 
throughout the calculation. 

In section y we discuss the shell-model Hamiltonian with a bound center-of-mass, the 
Faddeev equation, and the methods used to derive the starting-energy-independent effective 
interaction. Results of the calculations for the three-nucleon system are presented in section 
p| . In particular, we discuss the harmonic-oscillator frequency and the model-space-size 
dependences. Conclusions are given in section IIV. 



II. SHELL-MODEL HAMILTONIAN AND THE FADDEEV APPROACH 

In most shell-model studies the one- plus two-body Hamiltonian for the A-nucleon sys- 
tem, i.e., 

where m is the nucleon mass and V^ifi — fj) the nucleon-nucleon interaction, is modified 
by adding the center-of-mass harmonic-oscillator potential ^Amil'^R'^, R = ^ S^Li This 
potential does not influence intrinsic properties of the many-body system. It provides, 
however, a mean field felt by each nucleon and allows us to work with a convenient harmonic- 
oscillator basis. The modified Hamiltonian, depending on the harmonic-oscillator frequency 
Q, can be cast into the form 

A 
i=l 

The one-body term of the Hamiltonian can be re-written as a sum of the center-of-mass 

term i/^ = + |ylmn^^^, Pcm = SiLiPi; and a term depending on relative coordinates 
only. In the present application we use a basis, which explicitly separates center-of-mass 
and relative-coordinate wave functions. Therefore, the center-of-mass term contribution is 
trivial and will be omitted from now on. 
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For a three- nucleoli system, i.e., ^4 = 3, the following transformation of the coordinates 



Un-f2), (3a) 
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y=^i[Wi+r2)-r3], (3b) 
and, similarly, of the momenta 

P = - P-i) , (4a) 

q=\lli.Pi+ P2) - \/|p3 , (4b) 

can be introduced that brings the relative-coordinate part of the one-body harmonic- 
oscillator Hamiltonian into the form 

Eigenstates of this Hamiltonian, 

\nlsjt,XCSjT,JT) , (6) 

are then used as the basis for our calculation. Here n, I and A/", C are the harmonic-oscillator 
quantum numbers corresponding to the harmonic oscillators associated with the coordinates 
and momenta r,p and y,q, respectively. The quantum numbers s,t,j describe the spin, 
isospin and angular momentum of the relative-coordinate partial channel of particles 1 and 
2. S = ^ and t = | are the spin and isospin of the third particle, while is the angular 
momentum of the third particle relative to the center of mass of particles 1 and 2. The J 
and T are the total angular momentum and the total isospin, respectively. Note that for 
the ^H nucleus J = and T = ^. 

The Faddeev equation for the bound sytem can be written in the form 

H\4>) = , (7) 

with 

H = Ho + V{r)T . (8) 

Here, V{r) = 1/N(V2f) - |ml]2 

is the potential and T, which has the properties of a metric 

operator 0,^, is given by 

T = 1 + T(-) + , (9) 

with T^^^ and T*^~^ the cyclic and the anticyclic permutation operators, respectively. We 
derived a simple formula for the matrix elements of T^"-* + T^"*"-* in the basis (|^), namely 

{mhSiJitiMli^lSJlT, JT\T^-^ + T(+V2/2S2j2t2, ATs^a^JsT, JT) = -Sn,,N2 



X 



LS 




X 



;-i)^^+^^+*^+*^-'^^-'^(M/:ini/iL|n2/2Ar2/:2i^)d=3 + (m/iM/:i^|Ar2/:2n2/2^)d=3l , (10) 
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where iVj = 2nj + + 2J\fi + Ci,i = 1, 2, j = y2j~+T and {MiCiniliL\n2l2J^2^2L) d=3 is the 
general harmonic-oscillator bracket for two particles with mass ratio 3 as defined, e.g., in 
Ref . , where a compact formula is also given for calculating the brackets. The expression 
(|T0|) can be derived by examining the action of T^^^ and T^~^ on the basis states @. A 



similar derivation for a different basis is described, e.g., in Ref. JIT . 

Note that the eigensystem of the metric T (^) consists of two subspaces. The first 
subspace has eigenstates with the eigenvalue 3, which form totally antisymmetric physical 
states, while the second subspace has eigenstates with the eigenvalue 0, which form a not 
completely antisymmetric, unphysical subspace of states. It is possible to hermitize the 
Hamiltonian (|) on the physical subspace, where it is quasi- Hermit ian (see the discussion of 
quasi- Hermitian operators, e.g., in Ref. |P) . The Hermitized Hamiltonian takes the form 

H = Ho + T^/W{r)T'^'' , (11) 

where T operates on the physical subspace only. 

Apparently, the interaction V^r) is diagonal in the quantum numbers Af, C, J (and also 
in s,j, t due to the properties of the nucleon-nucleon potential). The metric T (^ is, on 
the other hand, diagonal in = 2n + / + 2A^ + L. Note that any basis truncation other 
than one of the type N < A^max violates, in general, the Pauli principle and mixes physical 
and unphysical states. Here, A^max characterizes the maximum of total allowed harmonic- 
oscillator quanta in the model space and is an input parameter of the calculation. At the 
same time, the truncation into totally allowed oscillator quanta < A^max preserves the 



equivalence of the Hamiltonians (|]) and (|TT|) on the physical subspace. 

From solving two-nucleon systems in a harmonic-oscillator well, interacting by soft-core 
potentials, one learns that excitations up to about SOOfiQ (A^max = 300) are required to 
get almost exact solutions. We anticipate, therefore, that at least the same number of 
excitations should be allowed to solve the three-nucleon system using the formalism discussed 
above. The Faddeev formulation has the obvious advantage compared with the traditional 
shell-model approach that the center-of-mass coordinate is explicitly removed. Even then, 
it is presently not feasible to solve the eigenvalue problem either for (H) or for (0) in 
such a large space. On the other hand, shell-model calculations are always performed by 
employing effective interactions tailored to a specific model space. In practice, these effective 
interactions can never be calculated exactly as, in general, for an A-nucleon system the A- 
body effective interaction is required. Consequently, large model spaces are desirable. In that 
case, the calculation should be less affected by any imprecision of the effective interaction. 
The same is true for the evaluation of any observable characterized by an operator. In the 
model space, renormalized effective operators are required. The larger the model space, the 
less renormalization is needed. We may take advantage of the present approach to perform 
shell-model calculations in significantly larger model spaces than are possible in conventional 
shell-model approach, particularly when using a Hermitized Hamiltonian (|11|). At the same 
time we can investigate convergence properties of effective interactions. If convergence is 
achieved, we should obtain the exact solution, since we recover the original full-space problem 
as A^max —>■ oo, provided that the condition Ves V is satisfied in this limit. 

Usually, the effective interaction is approximated by a two-body effective interaction de- 
termined from a two-nucleon system. In the present calculations we replace matrix elements 
of the potential V{r) by matrix elements of an effective two-body interaction, derived in 
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a straightforward manner for each relative-coordinate partial channel, 
nucleon Hamiltonian is then 



H2 = Hq2 + V 



-0 



2m 2 



The relevant two- 



(12) 



which can be solved as a differential equation or, alternatively, can be diagonalized in a 
sufficiently large harmonic oscillator basis. The latter possibility is, obviously, not applicable 
for hard-core nucleon-nucleon potentials. 

To construct the effective interaction we employ the Lee-Suzuki similarity transfor- 
mation method, which gives the effective interaction in the form PVesP = PVP + PVQuP, 
with uj the transformation operator satisfying u = QuP, and P, Q = 1 — P, the projectors 
on the model and the complementary spaces, respectively. Our calculations start with exact 
solutions of the Hamiltonian (0) and, consequently, we construct the operator u and, then, 
the effective interaction directly from these solutions. Let us denote the relative-coordinate 
two- nucleon harmonic-oscillator states, which form the model space, as \ap), and those 
which belong to the Q-space, as lag). Then the Q-space components of the eigenvector \k) 
of the Hamiltonian (0) can be expressed as a combination of the P-space components with 
the help of the operator u 



(13) 



ap 



If the dimension of the model space is dp, we may choose a set K. of dp eigenevectors, 
for which the relation ( P^ will be satisfied. Under the condition that the dp x dp matrix 
{ap\k) for I A;) G /C is invertible, the operator u can be determined from (0). In the present 
application we select the lowest states obtained in each channel. Their number is given by 
the number of basis states satisfying 2n + I < A^max- Once the operator u is determined the 
effective hamiltonian can be constructed as follows 



{-fp\H2es\ap) = 



{'^p\k)Ek{k\ap) + Y{-ip\k)Ek{k\aQ){aQ\uj\ap) 



(14) 



It should be noted that P\k) = J2ap \c(p){c(p\k) for \k) G /C is a right eigenvector of (|T3|) 
with the eigenvalue E^. 

This Hamiltonian, when diagonalized in a model-space basis, reproduces exactly the set 
/C of dp eigenvalues E^. Note that the effective Hamiltonian is, in general, quasi-Hermitian. 
It can be hermitized by a similarity transformation determined from the metric operator 
P(l + uj'^uj)P. The Hermitian Hamiltonian is then given by |T3[ 



2cS 



P{l + tu^u)P 



1/2 



2cS 



P(l + UJ^UJ)P 



-1/2 



(15) 



Finally, the two-body effective interaction used in the present calculations is determined 
from the two- nucleon effective Hamiltonian (O) as Ves = -f^2eff — -^^02- 



III. APPLICATION TO THE THREE-NUCLEON PROBLEM 

In this section we discuss the results of application of the formalism outlined in section 
H for the ^H system. In the calculations we use the Nijmegen II nucleon-nucleon potential 
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1^ corrected in the ^Pi wave [|T5|, and the Reid93 nucleon-nucleon potential |T^. We work 
in the isospin formahsm; the charge invariant potential Vn = |Kin + l^np is used for each 



T = 1 wave [|T§. 

The two-body effective interaction employed in the calculation is derived from the Eqs. 
(p!3|)-(p^). Our model space is characterized by the condition < A^max, N = 2n+l+2Af+C 
The condition for the relative-coordinate effective-interaction model space is then 2n + I < 
Nraa,x- When diagonalizing the two-nucleon relative-coordinate Hamiltonian (^) in the full 
space we truncate the harmonic-oscillator basis by keeping only the states with n < 152. The 
error caused by this truncation can be estimated, as the system can be solved as a differential 
equation. We found that the low-lying eigenvalues obtained in the two calculations do not 
differ by more than ^ 10^^ MeV and in most cases by much less. The lowest eigenvalues 
are typically of the order of 10^ MeV. Note that this error decreases with increasing Q. We 
calculated the effective interactions up to A^max = 32, as required in the present aplication. 

Once the effective interaction is found we may directly diagonalize the non-Hermitian 
Hamiltonian (|[) in the basis (§) truncated by A^max and with V replaced by V^g. 

On the other hand, a calculation with the hermitized Hamiltonian (|1T]) can be performed 
in three steps. First the effective interaction is calculated for each relative-coordinate partial 
channel. Second the metric T is diagonalized for each A^ up to A'max- The physical 
eigenvectors corresponding to the eigenvalue 3 are selected and used, finally, as a new basis 
in which the Hamiltonian ( pT]) is diagonalized. As the number of physical states is about 
a third of the number of all original basis states (j^), it is more efficient to diagonalize the 
Hamiltonian ([TT| ) than the non-Hermitian Hamiltonian (j^), in particular for higher values 
of A^max- In fact, for A^max > 22 we used the Hermitian Hamiltonian only. Note that the 
Hamiltonians (P) and (|lT]) have identical spectra of the physical states, provided that no 
other truncation than A^ < A^max is allowed. The unphysical eigenstates of H (j^) have 
energies corresponding to the unperturbed harmonic oscillator, starting at 3hQ. 

In Figs. ffl-H we present the results for the ground-state energies and point-nucleon radii. 



calculated from (r^) = — R)'^), obtained with the Nijmegen II and the Reid93 

nucleon-nucleon potentials, respectively. Our calculation starts at A'max = 8, which corre- 
sponds to a model space easily accessible with the traditional shell-model calculations. In 
Ref. 0] we performed an 8hQ calculation for ^H using a slightly different effective interaction 
than we employ here but derived in an analogous way. Note that it is straightforward to 
transform the relative-coordinate effective interaction used in the present calculations to the 
two-particle basis used for the shell-model input by the standard transformation [O. We 



used the transformed interaction in the 8hfl space to test our results. The shell-model di- 
agonalization was performed by using the Many-Fermion-Dynamics Shell-Model Code 



and we obtain the same answers from both the present calculation and the shell-model 
calculation. The present calculation has, obviously, much smaller dimension. 

As the results depend on A^max and fl introduced in Eq. (Q), we must test the convergence 
with regard to both of these parameters. With increasing A^max the calculations grow tedious. 
We performed the calculations up to A^max = 32 for a wide range of the harmonic-oscillator 
frequencies Q with values typical for standard shell-model calculations varying from hfl = 14 
MeV to hQ = 24 MeV. In general, we observe a slow convergence with increasing A'max- 
An unusual feature is the convergence from below. This is caused by the use of effective 
interactions instead of the free nucleon-nucleon interaction. The effective interactions we 



6 



employ are too strong. We have not reached the convergence with respect to Q in the 
whole range studied. However, for the values hfl = 22 — 24 MeV our results almost reach 
convergence with A^max = 32, in particular for the Reid93 potential. We note that the 
traditional 34-channel Faddeev calculation, as reported in Ref. gives the binding energies 
7.62 MeV and 7.63 MeV for the Nijmegen II and the Reid93 nucleon-nucleon potentials, 
respectively. We present these values in Figs. |l|, and |^ as dotted lines for comparison. From 
the figures it is apparent that we are obtaining virtually the same values in the calculations 
which start to converge. When comparing the results for the two different potentials used we 
can see larger sensitivity to Q of the Reid93 calculation. On the other hand, the calculation 
using the Nijmegen II potential is slower in reaching the stability with respect to N^^^^. 
While the ground-state energy calculation begins to stabilize for the largest values of fl 
employed, the radius calculation has not reached complete stability for any of the Q values 
within the model spaces we used. 

As a further test on the stability and convergence of the method we analyzed the ground- 
state wave functions and calculated the probability of S, P, D, S', S" states. In Fig. ^we show 
the D-state and S"-state probabilities as a function of A/'max for the Reid93 calculations with 
hfl = 19 and 24 MeV. We observe a good stability with respect to A^max and little dependence 
on Q for larger model spaces. The D-state probability approaches 8.4%, S"-state probability 
percentage reaches 1.2%. Not shown in the figure are the calculated P-state and ^"'-state 
percentage probabilities, for which we get 0.06%, and ~ 10~^%, respectively. The present 
numbers are in agreement with those obtained using other nucleon-nucleon potentials [Q. 
The D-state percentage is approximately 1.5 times the corresponding D-state percentage of 
deuteron (5.7%). 

In addition to the calculations discussed so far, we also computed properties of ^He, 
with the focus on obtaining the binding-energy difference between and ^He. In those 
calculations the Coulomb potential was added to the proton-proton potential and the av- 
eraged potential = f^p + |Kip was eventually used for each T = 1 wave [|TB|. The 
binding-energy differences obtained using the Reid93 potential and hQ = 19 and 24 MeV 
are presented in Fig. For larger model spaces we get an almost fi-independent difference. 
The binding-energy splitting shows convergence with increasing A'max- It decreaces with 
A^max in correlation with increasing point-nucleon radius and approaches 0.66 MeV. This 
result is again in agreement with those obtained using other nucleon-nucleon potentials . 
Note that the experimental value of the binding-energy difference is 0.764 MeV. To test 
the quality of the approximation used for the potential averaging and limitation to T = |, 
we performed an A^max = 8 calculation with complete isopin breaking using proton-neutron 
formalism. The present method can be used to perform calculations with isospin breaking. 
For this particular calculation, however, we employed the Many-Fermion-Dynamics Shell- 
Model Code. The effective interaction was calculated separately for the proton-proton and 
proton-neutron systems, respectively, and transformed to the two-particle basis as discussed 
earlier in this Section. In this way we found that the binding energy obtained with and 
without isospin breaking differs by 11 keV and the nucleon radius differs by less than 0.001 
fm in a calculation with hfl = 19 MeV. This confirms, that the limitation to T = | together 
with the potential averaging provides an excellent approximation. 

We stressed before that no other basis truncation than A^ < A^max was used. That 
means that we keep all the relative-coordinate channels in the basis. In most calculations, 
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however, we set the nucleon-nucleon potential to zero for j > 6. We also performed 
calculations with T^i set to zero for j > 4. The largest contribution of the j = 5,6 waves 
to the binding energy we observed was about 5 keV. This contribution increases with A^max 
and Q. Moreover, we performed several calculations with non-zero up to j = 9, and 
found that the ground-state energy is affected by less than 0.3 keV compared to the j < 6 
calculations. Note that such a truncation of does not imply that V = 0, see Eq. (|[). 
Also, this type of potential truncation is not the same as used in the traditional Faddeev 
calculations 0,^^,0. The difference is, that there is no truncation in the treatment of T 
(^) in the present calculations. 

It should be noted that the calculated binding energy obtained from calculations employ- 
ing nucleon-nucleon potentials fitted to the two-nucleon scattering underbind by about 
0.8 MeV, as its experimental binding energy is 8.48 MeV. Suggested solutions to this prob- 
lem include the use of three-body forces, non-local potentials and relativistic corrections 
0. 

IV. CONCLUSIONS 

In the present paper we have discussed the three-nucleon bound system solution by 
combining the shell-model approach with the Faddeev method. The use of Faddeev's de- 
composition reduces the basis and allows to perform shell-model calculations in significantly 
larger model spaces than in the traditional shell-model approach. We were able to calcu- 
late with the model spaces which included up to 32hQ (A^max = 32) harmonic-oscillator 
excitations. 

We employed effective interactions, which take into account the two-body correlations, in 
the calculations. These effective interactions were derived in the two-particle relative coordi- 
nate channels from the Nijmegen II and Reid93 nucleon-nucleon potentials and subsequently 
used in the three-body calculation. 

As our results depend on the model-space size parameter A'max and on the harmonic- 
oscillator frequency Q, we tested the convergence in both these parameters. Even for the 
largest model spaces we have not reached complete convergence with respect to Q in the 
whole range of the used values. However, for HQ = 22 — 24 MeV our results start to converge 
to the binding energies obtained in the standard Faddeev calculations. As we include more 
partial waves in the nucleon-nucleon potential, typically up to j = 6, our results seem 
to confirm the statements in Ref. |^ that the 34-channel standard Faddeev calculation 
converged within 0.01 MeV. We have seen, in fact, that the higher partial waves of the 
nucleon-nucleon interaction are not significant. However, we believe a proper treatment of 
the metric operator (j^) is important. 

We observed that the convergence was rather slow with some dependence on the type of 
the nucleon-nucleon potential. In smaller model spaces we cannot reproduce correctly the 
ground-state energy and the radius at the same time for any choice of Q. The wave function 
probability distribution shows good stability as well as the binding-energy difference between 
■^H and ^He. Our results show, that an 8hQ calculation accesible by the standard shell-model 
approach describes the ground states properties within 10% of exact values. Modifications 
of effective interactions are possible to improve the description in small spaces. Examples 
of such modifications can be found in Refs. I]^,^. 
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An unusual feature of the present approach is the convergence from below. This is 
caused by the use of effective interactions instead of the free nucleon-nucleon interactions. 
The effective interactions we employed were too strong. Obviously, it is possible to test 
convergence properties of alternative effective interactions, as discussed, e.g., in Ref. ||20|. It 



should be noted, however, that the effective interaction should meet the criterium Ves — V 
for A/'max — > C)0 in order to converge to the exact solutions. 

The formalism discussed here may be generalized for more complex systems as well. In 
particular, we are working on a generalization of the formalism for the A = 4 system, which 
relies on some results presented in this work. Also, it may be used to solve the three-nucleon 
system bound in a harmonic-oscillator well. Then, from those solutions three-body effective 
interactions can be constructed. Such interactions, after transformation to an appropriate 
three-particle basis can serve as an input to standard shell-model calculations for light nuclei. 
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FIGURES 



FIG. 1. Ground-state energy, in MeV, dependence on the maximal number of har- 
monic-oscillator excitation allowed in the model space for the Nijmegen II potential. Results 
for hfl = 14, 17, 19, 21, 22, 24 MeV are presented. The dotted line represents the result -7.62 MeV 
of the 34-channel Faddeev calculation reported in Ref. |||]. 

FIG. 2. Ground-state energy, in MeV, dependence on the maximal number of har- 
monic-oscillator excitation allowed in the model space for the Reid93 potential. Results for 
h^l = 14, 17, 19, 22, 23, 24 MeV are presented. The dotted hue represents the result -7.63 MeV 
of the 34-channel Faddeev calculation reported in Ref. H]. 

FIG. 3. Point-nucleon radius, in fm, dependence on the maximal number of har- 
monic-oscillator excitation allowed in the model space for the Nijmegen II potential. Results 
for nn = 14, 17, 19, 21, 22, 24 MeV are presented. 

FIG. 4. Point-nucleon radius, in fm, dependence on the maximal number of harmonic-oscillator 
excitation allowed in the model space for the Reid93 potential. Results for hil. = 14, 17, 19, 22, 23, 24 
MeV are presented. 

FIG. 5. D-state and S"-state probability, in %, dependence on the maximal number of har- 
monic-oscillator excitation allowed in the model space for the Reid93 potential. Results for 
h^l = 19, 24 MeV are presented. 

FIG. 6. Dependence of the energy difference, in MeV, between the binding energies of '^H and 
^He on the maximal number of harmonic-oscillator excitation allowed in the model space for the 
Reid93 potential. Results for HO, = 19, 24 MeV are presented. 
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